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ABSTRACT 

The possibility that long gamma-ray burst (GRB) jets are structured receives growing attention re- 
cently, and we have suggested that most GRBs and their softer, less energetic fraternity, X-ray flashes 
(XRFs), can be understood within a quasi- universal structured jet picture, given that the jet structure 
of each individual burst is a Gaussian or similar function. Here we perform a global test on such a 
quasi-universal Gaussian-like structured jet by comparing Monte-Carlo simulation results with a broad 
spectrum of observational data. Using the same set of input parameters as in the previous work (Zhang 
ct al. 2004), we confront the model with more observational constraints. These constraints include the 
burst redshift distribution, jet break angle distribution, two-dimensional redshift vs. jet break angle 
distribution, luminosity function, and logiV — logP distribution. The results indicate that the model is 
generally compatible with the data. This conclusion, together with our previous tests with the observed 
jet break angle vs. isotropic energy and observed peak energy vs. fluence relations, suggests that current 
long GRB and XRF data are generally consistent with such a quasi-standard-energy and quasi-standard- 
angle jet picture. With future homogeneous burst samples (such as the one to be retrieved from the 
Swift mission), the refined GRB jet structure can be further constrained through a global comparison 
between various observed and predicted burst property distributions and relations. 

Subject headings: gamma rays: bursts - jets 



1. INTRODUCTION 

The geometrical configuration is an essential ingredi- 
ent in characterizing and understanding astrophysical phe- 
nomena. There is growing evidence that long gamma- 
ray bursts (GRBs) are originated from collimated jets. 
This has been mainly suggested by an achromatic steep- 
ening break observed in many GRB afterglow light-curves 
(Rhoads 1999; Kulkarni et al. 1999; Harrison et al. 
1999). This interpretation receives indirect support from 
the intriguing fact that the geometry-corrected total en- 
ergy in the GRB fireball is essentially constant (Frail 
et al. 2001; Bloom, Frail, & Kulkarni 2003), i.e., E jet = 
Ei SO (l — cos9j) ~ const, where E iso is the total energy 
emitted in the gamma-ray band assuming isotropic emis- 
sion, and 9j is the jet angle inferred from the light-curve 
breaks. In view of these facts, there are two distinct ap- 
proaches in constructing jet models. One is that differ- 
ent GRBs collimate the same total energy into different 
angular openings, with the angular energy distribution 
within the jet being constant (Rhoads 1999; Frail et al. 
2001). Another is a family of quasi-universal structured 
jet models, where the structured jet has a power-law or 
a Gaussian (or functions in more general forms) angular 
energy distribution with respect to the jet axis (Zhang 
& Meszaros 2002a; Rossi, Lazzati, & Rees 2002; Lloyd- 
Ronning, Dai, & Zhang 2004). In the former scenario, the 
jet opening angle exclusively defines the jet break angle 6j , 
while in the later scenario, in most cases 0j is interpreted 
as the observer's viewing angle. It is essential to under- 
stand whether the GRB jets are structured, and if yes, how 
they are structured. This has important implications for 



the fundamental questions such as the total energy budget 
in the explosion, the physical origin of the collimation, and 
the birth rate of the GRB progenitor. 

In this paper, we define a "structured jet" as a jet with a 
certain functional angular distribution structure of energy 
(and possibly Lorentz factor as well), such as the power- 
law function with various indices, the Gaussian function, 
or numerous other possible structures one can think of. 
Since it is directly connected to the Ei SO oc 9J 2 (Frail et al. 
2001) correlation, one particular structured jet model, i.e., 
the power-law jets with index -2 (Rossi et al. 2002; Zhang 
& Meszaros 2002a), has received broad attention. We call 
this special type of jet as the "universal jet" following the 
convention in the literature since it has the potential to 
interpret all GRB data with a universal configuration. We 
notice that in some papers, "universal jets" and "struc- 
tured jets" have been used interchangeably which, to our 
opinion, may cause confusion to the readers, since a cer- 
tain criticism to the universal jet model may not apply to 
more general structured jet models with other jet struc- 
tures. The conventional top-hat jets are called "uniform 
jets" in this paper. 

The need of understanding GRB jets is boosted by the 
recent identification of X-ray flashes (XRFs, Heise 2003; 
Kippen et al. 2003), a fainter and softer version of GRBs, 
as a closely related phenomenon with GRBs. Recent obser- 
vations reveal another intriguing empirical correlation be- 
tween the cosmic rest frame GRB spectral peak energy and 
the isotropic gamma-ray energy, i.e., E pea k oc (Eiso) 1 / 2 , 
which was identified in the BeppoSAX (Amati et al. 2002) 
or even BATSE (Lloyd, Petrosian, & Mallozzi 2000) GRB 
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data, and was also found to extend to the XRF regime 
in the HETE-2 data (Lamb, Donaghy, & Graziani 2004; 
Sakamoto et al. 2004). Although obvious outliers (e.g. 
GRB 980425) exist, the correlation is found to be valid 
within the time-dependent spectra of invidual BATSE 
bursts (Liang, Dai, & Wu 2004a). This result strength- 
ens the empirical law, which suggests that it is related to 
some intrinsic physical processes. The relation is under- 
standable within the currently leading GRB models if a 
certain correlation between the bulk Lorentz factor and the 
burst luminosity is assumed (Zhang & Meszaros 2002b), 
i.e., r oc L k with k being different values for different mod- 
els. In particular, it is consistent with the internal shock 
model if the bulk Lorentz factor (and hence the internal 
shock radius) is insensitive to the burst luminosity, i.e., 
k = 0. In this paper, we assume that the E pea k — E iso 
correlation holds for the majority of GRBs and XRFs. 

In addition, the HETE-2 sample bursts also indicate 
another interesting fact that the contributions to the to- 
tal number of bursts from GRBs, X-ray rich GRBs (XR- 
GRBs), and XRFs are approximately equal (Lamb et al. 
2004). This fact presents an important criterion to test 
the validity of any jet model. 

The similarities between GRBs and XRFs have stim- 
ulated studies towards unifying the GRB and XRF phe- 
nomena through different geometrical configurations (e.g. 
Lamb et al. 2004; Yamazaki, Ioka, & Nakamura 2003; 
Zhang et al. 2004). Lamb et al. (2004) pointed out 
that the universal jet model is inconsistent with the GRB- 
XRGRB-XRF number ratios detected by HETE-2, and 
turned to suggest a uniform jet model for all GRBs and 
XRFs. If one assumes a standard energy budget for all 
GRBs and XRFs, such a uniform jet model unavoidably 
leads to the conclusion that GRBs have very narrow jets 
with typical opening angle smaller than 1 degree. The 
universal jet model was also tested against various criteria 
recently (e.g. Perna, Sari, & Frail 2003; Nakar, Granot, 
& Guetta 2004; Guctta, Piran, & Waxman 2004a), and it 
has been found that such a model may violate some obser- 
vational constraints. However, a pure universal jet model 
corresponds to a strict Ei so cx 9J 2 relation. In reality, the 
data indicate that this correlation is only valid in a statis- 
tical sense (Frail et al. 2001; Bloom et al. 2003). In the 
Ei S o — 6j plane, the afterglow data are distributed around 
the E iso cx 9~ 2 line with moderate scatter (Lloyd-Ronning 
et al. 2004). This fact alone already suggests that GRB 
jets are not "universal" . Any jet model aiming to interpret 
the GRB phenomenology is at best "quasi- universal" , i.e., 
different jets may share a more or less the same structure, 
but the parameters to define the structure should have 
some scatter around some typical values. An important in- 
sight is that when parameter scatter is taken into account, 
the jet structure is no longer obliged to be power-law with 
a -2 index. Other jet structures are also allowed (Lloyd- 
Ronning et al. 2004). In particular, we (Zhang et al. 2004) 
recently proposed a quasi-universal model for GRBs and 
XRFs. In order to successfully reproduce the right rela- 
tive numbers of GRBs, XRGRBs, and XRFs, we suggest 
that the jet structure in individual bursts is Gaussian-like, 
or with a similar structure. This ansat was verified with a 
Monte-Carlo simulation, and the model can also reproduce 
the E iso — 6j relation. In view that the narrow uniform 



jet model (Lamb et al. 2004) conflicts with the standard 
afterglow model (Zhang et al. 2004), and that the uni- 
versal jet model encounters various difficulties (e.g. Lamb 
et al. 2004; Guetta et al. 2004a), we tentatively suggest 
that the quasi-standard-cncrgy and quasi-standard-anglc 
Gaussian-like jet model is a more plausible one to inter- 
pret GRB and XRF data in a unified manner. 

In order to prove this suggestion, the quasi-universal jet 
model needs to confront a broader spectrum of data. Since 
within a structured jet model the probability of observing 
the jet at angle 9 V is proportional to sin^), many ob- 
servational properties can be predicted once the jet struc- 
ture function and the variation parameters are given. In 
this paper, besides the E° h e s ak vs. fluence relation and the 
Ei 8 o — 0j relation we have already tested in Zhang et al. 
(2004), we consider several new constraints including the 
burst redshift (z) distribution, jet angle (9j) distribution, 
two-dimensional z — 9j distribution, luminosity function, 
and log N — log P distribution. Some of these criteria have 
been taken individually to test some jet models (e.g. Perna 
et al. 2003; Lloyd-Ronning et al. 2004; Lin, Zhang, & Li 
2004; Liang, Wu, & Dai 2004b; Nakar et al. 2004; Guetta 
et al. 2004a). However, none of the previous studies per- 
formed a global test for a particular model with all the 
criteria. We believe that such a global test is essential to 
constrain and to finally pin down the right GRB jet struc- 
ture. Here we perform such a test with the quasi-universal 
Gaussian- like jet model (Zhang et al. 2004). The exact 
GRB structure may differ from the simple Gaussian form. 
We take this simple structure as the starting point to ex- 
amine how well it could reproduce the data. Due to the 
"quasi-universal" nature, analytical studies may not be ad- 
equate, and we perform a set of Monte-Carlo simulations 
to access the problem. 

2. MONTE-CARLO SIMULATIONS 

We perform Monte-Carlo simulations for a quasi- 
universal Gaussian-like jet model. The jet structure and 
the input parameters that we use are the same as those 
used in Zhang et al. (2004), where the motivation to in- 
troduce such a jet structure is also explained. Below we 
describe the parameters of this model in more detail. 

First, we approximate the angular energy distribution 
of the jet as 

_ e 2 

e{9) = e e (I) 

The total energy of the jet, Ej, is obtained by integrating 
e{9) over the entire solid angle 

,tt/2 

Ej = 4tt / e{9) sin(9)d9. (2) 
Jo 

With a small 6 , the total jet energy is approximately 
Ej ~ 2ne 9 2 (Zhang & Meszaros 2002a). This jet struc- 
ture contains two parameters, the total energy of the jet, 
Ej, and the characteristic jet width, 9$. The parameters, 
Ej and 6>o, are distributed in log-normal distributions for 
the simulated bursts. This quasi-universal approach is re- 
quired to reproduce the large scatter of the E iso -9j relation 
of GRBs (Lloyd-Ronning et al. 2004). In particular, these 
two parameters are constrained to be around (Zhang et al. 
2004) 

<log(JL)> „ 51.1 (3) 
I erg 
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Although the definition Ej in the Gaussian structured jet 
model (Eq. [2]) is different from that in the uniform jet 
model, our best fit typical jet energy (Eq. [3]) is consis- 
tent with the one in the uniform jet model (Bloom et al. 
2003). We have shown that this set of input parameters 
can roughly reproduce the approximately equal numbers 
of GRBs, XRGRBs, and XRFs, the E iso -9j relation, and 
the -E°g* fe -fluence relation. 

Second, in a structured jet model the observing angle, 
0,,, is distributed as 



dN(0 v ) 
d6 v 



= sin(0„). 



(7) 



The isotropic equivalent energy, E iso , is defined as 

E tso = 4ire(9 v ). 

The jet break angle, 6j, of a Gaussian jet is (Kumar & 
Granot 2003; Zhang et al. 2004; Rossi et al. 2004) 
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Finally, the number of bursts per unit redshift, N(z), is 
distributed as 



dN(z) _ R GRB (z)dV(z) 



(10) 



dz 1 + z dz 

where dV(z)/dz is the comoving volume per unit redshift 
and Rgrb is the GRB rate. The comoving volume is ob- 
tained as 

4 ^ C [o ro (i + zf + n k (i + zf + ^a]- 1 / 2 , 

(ii) 

where Ho = 70 km s 1 Mpc 1 , Cl m — 0.3, fife = 0, and 
Oa = 0.7, and Dl is the luminosity distance. We assume 
that the GRB rate traces the star forming rate such that 
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Here, we used the Rowan-Robinson star forming rate 
(Rowan-Robinson 1999, cf. Lin et al. 2004). Two sets of 
burst redshifts are simulated with the parameter z peak = 2 
and z pea k — 1, respectively. 

We simulate 10,000 bursts. For each burst we simu- 
late Ej, 9q, 9 v , and z according to the distributions those 
parameters follow. We calculate other parameters of the 
simulated bursts in the following equations. The peak en- 
ergy of the simulated bursts are calculated through the 
E iso -E pea k relation (Amati et al. 2002) 



E pe ak ~ 100 keV( 



Eiso 

10 52 erg 



1/2 



(13) 



in which we introduced lognormal scatter (with a ~ 0.3) 
to reflect the statistical nature of the observed data points. 
The observed peak energy is related to the rest-frame peak 
energy by 

E; b eak = E peak /(l + z), (14) 



and the burst energy fluence is calculated as 

E iso (l + z) 



F = 



(15) 



In order to obtain the peak flux or peak luminosity of the 
bursts, we simulate the conversion time scale T in the rest 
frame of the bursts as a log-normal distribution (see also 
Lamb et al. 2004), 



T 



- 0.56 
0.05 , 



such that 
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fpeak 



Ei S o/T 
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(16) 
(17) 

(18) 
(19) 



T(l + z) 

The central value of this distribution, 3.63 s, is consistent 
with the value of 3.41 s obtained from the HETE-2 and 
BeppoSAX burst sample (Lamb et al. 2004). We note 
that this conversion time scale should be shorter than the 
true duration of the burst since the peak flux is higher 
than average flux. The exact central value used in this 
paper is obtained as the best fit to the BATSE peak flux 
distribution through a series of simulations with different 
central values (please see § 3.5). 

The simulated bursts are filtered by different detection 
thresholds that represent the sensitivities of different sur- 
veys. In principle, the thresholds should be calculated in 
unit of peak photon flux. The peak photon flux depends 
on the spectral shape of the bursts, and the sensitivity is 
also a function of photon energy. Moreover, most of the 
bursts observed are strongly variable. For simplicity, here 
we use the simulated burst fluence to define various detec- 
tion thresholds in most of the simulations except for the 
log N — log P distribution. 

3. RESULTS 

After obtaining the sample of the simulated bursts with 
various burst parameters calculated, we compare the sim- 
ulation results with current observations and perform a 
global test to the GRB jet structure with the observational 
constraints. The constraints include the burst redshift 
distribution, jet break angle distribution, redshift vs. jet 
break angle two-dimensional distribution, luminosity func- 
tion, log A" — logP distribution, E iso - 9j relation, E° b e s ak 
- fluence relation, and the relative numbers of GRBs, XR- 
GRBs, and XRFs. We note that the last three tests have 
been presented in a previous paper (Zhang et al. 2004). 
Here, we only present the rest of the simulation results. 

3.1. Redshift Distribution 

There are two major uncertainties for the simulated 
burst redshifts. First, we assume that the GRB rate trace 
the star forming rate which may involve some uncertain- 
ties. In addition, the underlying star forming rate itself is 
not accurately constrained. The observed redshift distri- 
bution of the bursts can be used, as a first step, to con- 
strain the star forming rate used in the simulation. We 
note that the jet structure will also affect the shape of the 
redshift distribution due to the limited sensitivities of the 
detectors. A homogeneous sample would be most suitable 
for these studies. The current sample of the bursts with 
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rcdshift measurements is, however, small and inhomoge- 
neous. Nonetheless, we compare this sample with the sim- 
ulation results to obtain some preliminary understanding 
of the issue. 

We plot the rcdshift distributions of the simulated bursts 
and compare with the sample obtained from the obser- 
vations (Figure 1). Figures la and lb represent simula- 
tions with two different star forming rate, with z pea k = 2 
and Zp ea k = 1) respectively. The observed burst redshifts 
are obtained from Bloom et al. (2003) and the GRB Lo- 
calization website maintained by J. Grciner, and there 
are 37 bursts with measured redshifts in total. We se- 
lect the simulated bursts with different detecting thresh- 
olds. The thresholds are set to be l.OxlO -5 erg cm -2 and 
5xl0 -8 erg cm -2 . The first threshold is selected to be 
much higher than the triggering thresholds for most de- 
tectors in order to reflect the selection effect that only a 
small number of bright GRBs have their redshifts mea- 
sured. For example, the median fluence value of the 28 
bursts with redshift measurements in Bloom et al. (2003) 
is 2.3xl0 -5 erg cm -2 . The second threshold is selected to 
match the sensitivity of HETE-2. 

Figure 1 shows that the observed redshift distribution 
peaks at z ~ 1. The observed distribution is consistent 
with both simulated distributions with different z pea k, if a 
high detection threshold of l.Ox 10 -5 erg cm -2 is adopted. 
This high threshold is reasonable to apply here when com- 
paring the simulated redshift distribution with the cur- 
rently observed redshift distribution because only a few 
bursts have their redshifts measured. These bursts only 
account for a small fraction of the total amount of bursts 
detected, and they are typically brighter. In order to bet- 
ter illustrate this point, we show in Figure 2 the redshift 
vs. fluence plots for the simulated bursts and the ob- 
served bursts from Bloom et al. (2003) sample. Both 
star forming rate models have been plotted. We can see 
that most of the observed bursts are concentrated in the 
high fluence end compared with the simulated bursts. The 
current observed redshift sample cannot distinguish be- 
tween the two star forming rates. However, as a more 
homogeneous GRB redshift sample is accumulated (with 
future instruments such as Swift), the peak of the GRB 
rate can be constrained as the two models predict differ- 
ent peaks as the detection threshold decreases. We adopt 
Zpeak = 2 when discussing the simulation results below. 
An even larger and more homogeneous sample is needed 
in order to constrain the shape of star forming rate, es- 
pecially whether the star forming rate is above or below 
the Rowan- Robinson rate at large redshifts (z > 2). For 
Zpeak = 2, the current data indicate tentative evidence 
that the nearby events are more abundant than what is 
expected from the model that assumes the standard star- 
forming rate. For z pea k — 1, the current model meets the 
data in the low-redshift regime as well. Whether there 
exist extra nearby GRBs is of great interest in the GRB 
community, and more redshift data are needed before a 
firm conclusion is drawn. 

We note that most XRFs do not have redshift measure- 
ments, which may present a bias in the observed sample. 
However, this should not affect the result too much be- 
cause the current redshift sample of GRBs is small (37 

4 The website is at http://www.mpe.mpg.de/~jcg/grbgcn.html 



in total) . Considering the number ratio of GRBs to XRFs 
from HETE-2 data, about 12 XRFs should be added to the 
sample. If the rcdshift distributions of GRBs and XRFs 
are similar, the addition of about 12 XRFs should not 
change the shape of the redshift distribution very much. In 
addition, the best fits to the observed redshift distribution 
are the simulated bursts with large fluences, and accord- 
ing to the -E°g* fc -fluence relation of the HETE-2 bursts, 
these high fluence bursts (> l.OxlO -5 erg cm -2 ) should 
all be GRBs. Similar arguments also apply to the follow- 
ing simulations of the jet break angle distribution and the 
two-dimensional jet break angle vs. redshift distribution. 

3.2. Jet Break Angle Distribution 

We plot the distribution of the jet break angle of the 
simulated bursts with different detection thresholds and 
the observed jet break angle distribution in Figure 3. The 
observed jet break angles are obtained from Bloom et al. 
(2003), and there are 16 bursts with jet break angle mea- 
surements. We exclude the bursts with upper or lower 
limit measurements on the jet break angles when compar- 
ing with the simulation results. The shape of the observed 
jet break angle distribution is sensitive to the bin size cho- 
sen because of the small sample size. However, the peak 
of the jet angle distribution can roughly be constrained at 
about 7 degrees. In this simulation, we also apply different 
fluence filters to the simulated bursts to simulate the ef- 
fect of the limited detector sensitivities. In particular, we 
adopt the fluence thresholds of l.OxlO -5 , 5.0xl0 -7 , and 
5.0xl0 -8 erg cm -2 , respectively. Figure 3 shows that the 
simulated jet break angle distribution is also sensitive to 
the detection threshold used in the simulation. The peak 
of the simulated distribution will move to larger angles 
when a lower threshold is selected. This is consistent with 
the result of Perna et al. (2003), who discovered this effect 
from the universal jet model. The simulated distribution 
with a fluence limit of l.OxlO -5 erg cm -2 is consistent 
with the observed jet break angle distribution. Consid- 
ering the difficulties in identifying the optical afterglows 
and in measuring the jet break angles, such a high thresh- 
old is a reasonable choice when comparing the simulation 
with the observations. The median fluence of the 16 bursts 
with jet break angle measurements is 2.3xl0 -5 erg cm -2 
(Bloom ct al. 2003), even higher than the highest fluence 
limit we adopt. Only five of the 16 bursts have fluences 
lower than l.OxlO -5 erg cm -2 , and the faintest one has 
a fluence of 3.17xl0 -6 erg cm -2 . Considering that there 
are some bursts that are brighter than the fluence limit we 
adopt but whose jet break angles are still not measured, it 
is more reasonable to compare the median fluence, rather 
than the smallest fluence, of the 16 bursts. In Figure 4, we 
show the fluence vs. jet break angle plot for the simulated 
bursts and the observed bursts (Bloom et al. 2003). Most 
of the observed bursts are in high fluence regions. We also 
compare the jet break angle distribution with the predic- 
tions from the z pea k = 1 Rowan-Robinson star forming 
rate, and lead to similar results. 

We reach a similar conclusion as Perna et al. (2003) that 
the predicted jet break angle distribution of structured jets 
is consistent with the currently observed sample distribu- 
tion. While Perna et al. (2003) discussed the universal 
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jets, our simulations are for Gaussian-like jets. We note 
that the detection threshold used in Perna et al. (2003) is 
the 90% efficiency peak flux threshold for BATSE, which is 
much more sensitive than the threshold we have adopted in 
the simulation that results in consistency between the data 
and the simulation. Since only a small fraction of bright 
GRBs have jet break angle measurements, it may be more 
appropriate to use a higher threshold than the BATSE 
detection threshold. When a higher threshold is selected, 
the 8j distribution peak for the universal jet model should 
move to a smaller value comparing with the observed one. 
Liang et al. (2004b) also noticed this independently and 
also adopted a high threshold in their simulations. 

3.3. Redshift vs. Jet Break Angle 

The jet break angle distribution discussed previously is 
a one-dimensional distribution which includes bursts at 
all redshifts. As pointed out by Nakar et al. (2004), a 
more accurate test is to perform a two-dimensional (z — 9j) 
distribution comparison between the data and the model 
prediction. It is possible that the two-dimensional dis- 
tribution does not agree with the observations while by 
integrating over redshift the onc-dimcnsional distribution 
agrees with the observations by chance. 

In Figure 5, we plot the simulated data points with dif- 
ferent fluence thresholds (same as those used in the one 
dimensional plot) together with the 16 bursts with both 
redshift and jet break angle measured (Bloom et al. 2003). 
The density of the data points represents the probability 
density function (PDF) of this two-dimensional distribu- 
tion. The plot shows that the PDF depends on the thresh- 
old of the detector. As the threshold goes higher, the peak 
of the PDF moves towards the region containing smaller 
jet break angles. This is consistent with the result from the 
one-dimensional analysis. In particular, the distribution of 
simulated bursts with the highest threshold is consistent 
with the observational distribution. Even we limit the data 
points from a narrow redshift bin (0.8 < z < 1.7) that in- 
cludes most of the observational data points, the simulated 
distribution and observational distribution is still consis- 
tent. Again, a high threshold is reasonable in this analysis 
since the bursts with jet break angle measurements are 
much brighter on average than the ones in the whole sam- 
ple of the detected bursts. 

Nakar et al. (2004) argued that this two-dimensional 
distribution from the prediction of a universal jet does not 
agree with the observation, especially for the bursts within 
the redshift range of 0.8 < z < 1.7. As we have shown 
previously (Lloyd- Ronning et al. 2004), the E iso — 8j data 
already require that the model has to be quasi-universal, 
thus the inconsistency suggested by Nakar et al. (2004) 
is largely due to their adopting a non-varying universal 
jet model. In particular, the sharp boundary of the re- 
gion allowed by the universal model should smooth out 
if a quasi-universal model is considered (as shown in Fig- 
ure 5). Moreover, the detection threshold used in Nakar 
et al. (2004) is also the BASTE detection threshold, which 
is much more sensitive compared with the burst sample 
with jet break angle measurements. Also we adopted a 
Gaussian-like jet structure, while they stick to the power- 
law structure which suffers other problems (e.g. numbers 
of XRFs respect to GRBs and the log N — \ogP distribu- 



tion) as well. 

3.4. Luminosity Function 

We plot the luminosity function of the present quasi- 
universal Gaussian-like jet model in Figure 6. This lu- 
minosity function is an update of the result presented in 
Lloyd-Ronning et al. (2004). In Lloyd-Ronning ct al. 
(2004), the total energy of the jet was taken as a constant 
rather than quasi-universal. When the total energy scat- 
ter is introduced, the simulated luminosity function does 
not show a bump at the break of the powerlaw index (cf. 
Figures 9 and 10 in Lloyd-Ronning ct al. 2004). In addi- 
tion, we have performed a rigorous calculation of the total 
energy instead of using the approximation for small an- 
gles. At large angles, this difference is more than a factor 
of two. 

The simulated luminosity function can be characterized 
by a broken power law, with a power law index of ~ — 2 
at high luminosity end (L > 10 52 erg s _1 ) and an index of 
~ —1 for low luminosities (L < 10 52 erg s _1 ). The simu- 
lated luminosity function is consistent with previous simu- 
lations of Gaussian jets performed in Lloyd-Ronning et al. 
(2004), except that there is no big bump at the break of 
the power index this time. Currently, the GRB luminosity 
function is not directly determined from the observations, 
since the sample of bursts with redshift measurement is too 
small. Nonetheless, there are several attempts to constrain 
the luminosity function through various approaches (e.g. 
Schmidt 2001; Norris et al. 2002; Lloyd-Ronning, Fryer, 
& Ramirez-Ruiz 2002; Stern, Tikhomirova, & Svensson 
2002; Firmani ct al. 2004). In general, many of these 
studies found a break of the luminosity function, with the 
power index steeper in the high luminosity range and flat- 
ter in the low luminosity range (see detailed discussion in 
Lloyd-Ronning ct al. 2004 and references therein). This 
is consistent with the simulated luminosity function. In 
particular, the shape of the simulated luminosity func- 
tion in this paper is very similar with that obtained from 
Schmidt (2001), except for the low luminosity end below 
10 50 erg s _1 where the luminosity functiom from Schmidt 
(2001) turns over. This luminosity function obtained from 
Schmidt (2001) can fit the BASTE log A - logP quite 
well assuming a certain star forming rate (Schmidt 2003). 

3.5. log N — log P Distribution 

We use the log A — logP plot as a final test to the jet 
structure used in this paper. The log A — logP plot has 
been used to constrain the jet opening angle distribution 
for the uniform jet and the star forming rate in previ- 
ous studies (Lin et al. 2004; Guetta et al. 2004a). In 
particular, Guetta et al. (2004a) pointed out that the 
log A — log P distribution predicted by the universal jet 
model over predicts bursts with faint flux. This incon- 
sistency is another manifestation of the problem of over- 
producing XRFs in the universal jet model (Lamb et al. 
2004). Since the quasi-universal Gaussian jet model can 
overcome the latter difficulty, it is natural to expect that 
it can solve the former problem as well. 

We use the bursts from the offline re-analyzed BATSE 
catalog (Kommers et al. 2000) including both triggered 
and untriggered bursts. The catalog includes a total of 
2167 bursts, in which 1393 are triggered bursts and 874 
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are untriggered bursts. The bursts selected in this cata- 
log are all long GRBs, which are directly related to our 
model. In order to estimate the simulated peak photon 
flux, we assume a Band function (Band et al. 1993) for 
the simulated bursts, and adopt the low- and high-energy 
photon indices as a = — 1 and (3 = — 2. As our simula- 
tion stems from the data provided by Bloom et al. (2003) 
where the isotropic energy is given in the 20-2000 keV 
band rest frame, we need to take into account the differ- 
ence between the bandpass of BATSE and that used in 
Bloom et al. (2003). The fluence that we obtained from 
Eq.(15) should be 

/•2000/(l+z)koV 

F= EN(E)dE . (20) 

J20/(l+z)koV 

The photon fluence in the BATSE band (50-300 keV) 
reads 



pph 
BATSE 



300koV 



N(E)dE . 



(21) 



50koV 



We can then calculate the photon fluence in the BATSE 
band using the energy fluence, i.e. 

Fj^N(E)dE 



jpph 
BATSE 



(22) 



r 2000/(l+z)koV . , ' 

J 2 o/(i+z)kov bN{E)dE 

The ratio between the two integrals in Eq.(22) only de- 
pends on the value of E p (which is simulated in the code) 
when both a and (3 are assigned to their typical values 5 . 
Essentially, we use the energy fluence to determine the 
normalization of the Band function and calculate the pho- 
ton fluence with the determined Band function. Finally, 
the peak photon flux can be obtained by dividing P^tse 
by the conversion timescale T defined in Eq.(16), i.e., 

(23) 



P, 



rph _ * BATSE 

/BATSE T (l + Z) 



We do not perform a fluence or energy truncation for the 
log A — logP analysis (unlike in the previous sections). 
The BATSE sample is much more homogeneous, and is 
free of most of the selection effects encountered for the 
other samples we have discussed (the sample with red- 
shift or jet angle information). The simulated sample 
is natually truncated with the offline detection thresh- 
old in peak photon flux (Kommers et al. 2000), i.e., 
0.18 photon cm" 2 s _1 . This is lower than the BATSE 
onboard detection threshold, 0.3 photon cm~ 2 s _1 . This 
peak photon flux truncation is adopted when we compare 
the simulated log A — log P distribution with the observa- 
tion through the Kolmogorov-Smirnov (K-S) test. Guided 
by the K-S test, we fit the simulated log A — logP distri- 
bution to the observed distribution with a normalization 
parameter so that the simulated log A — log P plot can 
be shifted vertically. The normalization reflects the differ- 
ence between the number of the simulated burst and the 
number of the true bursts in the BATSE sample. After 
the fitting, the simulated bursts fit very well with the ob- 
served log N — log P plot. The log A — log P distributions 
for the simulated bursts after the fitting and the bursts 
from BATSE catalog are shown in Figure 7. We perform 
a K-S test to the two cumulative distributions after the 
peak photon flux truncation at 0.18 photon cm~ 2 s _1 and 

5 We have also tested the cases when a and deviate from the nomical values. The results essentially remain unchanged. 



a K-S chance probability of 11% is obtained, which impli- 
cates that the simulated log N — log P distribution cannot 
be rejected through the K-S test. 

We also test the model with the z pea k = 1 Rowan- 
Robinson star forming rate and similar results are ob- 
tained. This result suggests that the observed BATSE 
log A — logP distribution can be reproduced with this 
quasi-universal Gaussian structured jet model. We also 
compare the BATSE catalog from Kommers et al. (2000) 
with the GUSBAD catalog from Schmidt (2004), and the 
log A — log P distributions obtained from the two catalog 
are generally consistent. 

The simulated log A— log P plot can generally reproduce 
the turnover behavior at the faint end, but still deviates 
the BATSE data at the lowest peak photon flux of the dis- 
tribution. This does not affect the main result because this 
deviation occurs at regions where the peak photon flux is 
below the limit of the BATSE offline search sensitivity so 
that a significant part of the burst could be missed. In 
addition, at low flux region the distribution is sensitive to 
the detection threshold adopted. We plot the simulated 
distribution extending down to a lower peak photon flux 
limit than BATSE as a prediction for future observations. 
In the simulation we did not perform any luminosity or 
Epeak truncation so that XRFs are also included in the 
simulated samples for this test. It is unclear if most XRFs 
have been detected with BATSE. Since a Band function 
(Band et al. 1993) have a typical upper photon spectral 
index /3 ~ — 2 for GRBs, XRFs should be detected as faint 
GRBs as well if their high energy photon spectral indices 
have a similar value. In addition, we also test the case by 
excluding the simulated burst with lower observed peak 
energies (e.g. those with E° b e s ak < 10 keV), and the re- 
sulting log A — log P plot is still consistent with the data. 
The reason is that XRFs mainly contribute to the faintest 
population of the distribution. 

We note that the simulated log A — logP plot is also 
sightly different from the observed distribution at high 
photon fluxes. However, the difference is not significant 
because the number difference between the simulated and 
observed bursts is small. The difference could result from 
small number statistical effects. This could be tested by 
increasing the number of simulated bursts. It may also in- 
dicate that the jet structure is more close to a power-law 
structure at small angles. 

Guetta et al. (2004a) fitted the predicted log A - log P 
distribution from the universal jet model and compared 
it with the observations. The total bursts (595) used in 
Guetta et al. (2004a) is smaller than the burst sample 
(2167 in total) used in this paper. However, even with a 
small sample, Guetta et al. (2004a) concluded that the 
universal jet model can be rejected from the log A — logP 
distribution. With a quasi-universal Gaussian-like jet 
model, the observed log A— log P distribution can be fitted 
quite well from the simulations performed in this paper, es- 
pecially in the range of —0.7 < logP (photon cm~ 2 s _1 ) < 
1 where the log A — log P distribution is constrained most 
accurately. The major difference between a power-law jet 
and a Gaussian jet is at large angles. The exponential drop 
of the jet energy at large angles in the Gaussian model is 
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the key to reduce the right number of faint bursts in the 
BATSE sample. The problem faced by the universal jet 
model in the log N — log P distribution test is avoided with 
the quasi-universal Gaussian jet model. 

Recently Guetta, Granot, & Begelman (2004b) claims 
that a power law jet truncated at large angles can also 
reproduce the BATSE log A - logP distribution. This 
is consistent with our argument raised in this paper. In 
view that a Gaussian jet invokes a natural exponential 
drop off at large angles while a reasonable power-law jet 
model has to invoke both a small angle and a large angle 
artificial break, we deem that the quasi-universal Gaus- 
sian jet model is a more elegant one. We notice that 
the log A — logP distribution is more sensitive to the jet 
structure than to the GRB distribution with rcdshift. It 
is therefore a powerful tool to pin down the correct jet 
structure. 

4. CONCLUSION AND DISCUSSION 

We perform Monte-Carlo simulations of a quasi- 
universal Gaussian-like structured jet and compare the 
simulation results with a wide spectrum of current ob- 
servations. The simulation results for the Gaussian-like 
jet used in the paper are generally consistent with vari- 
ous observational constraints, including the burst rcdshift 
distribution, jet break angle distribution, two-dimensional 
distribution of rcdshift and jet break angle, luminosity 
function, and log A — logP distribution. This result is 
complementary to the previous simulation results (Zhang 
et al. 2004) that showed that the number ratios among 
GRBs, XRGRBs and XRFs, the observed jet break angle 
vs. isotropic energy relation, and the observed peak en- 
ergy vs. fluence relation are consistent with predictions 
from this jet model. 

Although the samples of different burst properties used 
for some tests in this paper are small so that very de- 
tailed constraints on the jet structure cannot be achieved, 
the global test performed on the quasi-universal Gaussian 
jet model (including this paper and the previous work, 
Zhang et al. 2004) at least suggests the following two ma- 
jor conclusions. First, in order to unify GRBs and XRFs 
through viewing angle effects such that GRBs are viewed 
from small observing angles and XRFs are viewed from 
large observing angles, the jet energy at large angles must 
drop exponentially. This comes from two constraints, i.e., 
the log A— log P distribution and the number ratios among 
GRBs, XRGRBs and XRFs. This requirement is con- 
sistent with the Gaussian jet model. The power-law jet 
structure cannot extend to large angles. Otherwise, it will 
over produce the number of XRFs and low flux bursts in 
the log A — log P distribution. Second, "quasi-universal" 
should also be an essential ingredient for the jet models, 



as it is unrealistic to assume that all the GRB progenitors, 
their environments, and other properties are exactly the 
same. The quasi-universal nature is crucial to fit the E iso - 
9j relation, and it could produce the luminosity function 
power-law index break that was implicated in many stud- 
ies. A quasi-universal Gaussian-like jet is suitable for both 
constraints. We also note that many burst property dis- 
tributions are sensitive to the detection threshold selected, 
and a suitable threshold for the current sample should be 
selected carefully when comparing the model predictions 
with the observations. 

Currently, the detailed jet structure at small angles is 
not well constrained. The slope may be steeper than that 
in the Gaussian model. It is possible that the Gaussian jet 
slightly under-predicts bright bursts as indicated from the 
^peofc _nuence Pl°t an d the E iso — 9j plot in Zhang et al. 
(2004). Such a deficit, however, could be well due to selec- 
tion effects since the brightest bursts are those most easily 
to detect and to localize. If the deficit is real, it may be 
accounted for through possible evolution of the GRB lu- 
minosity function, which was suggested by recent studies 
(Lloyd-Ronning et al. 2002; Wei & Gao 2003; Yonetoku 
et al. 2004; Graziani et al. 2004). It may also be under- 
stood in terms of the two-component jet picture (Berger 
etal. 2003; Huang et al. 2004; Liang & Dai 2004), where 
a bright core component contribute to more bright bursts. 
These issues can be tested thoroughly with a larger, ho- 
mogeneous sample of bursts accumulated. 

In the Swift era, the sample of GRBs with redshift and 
jet break angle measurements is anticipated to be much 
larger than the current sample. In addition, these bursts 
will form a homogeneous sample that is most suitable to 
apply statistical analyses on various burst properties. We 
anticipate most of the GRB relations used in this paper 
and Zhang et al. (2004) will be constrained more accu- 
rately, except for the log A — logP distribution. With 
refined Monte-Carlo simulations of this and other models, 
the details of GRB jet structure can be pinned down more 
precisely through statistical analyses of the observed burst 
properties and the model predictions. 
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Fig. 1. — (a) The probability distribution functions (PDF) of GRBs with respect to redshift for the observed bursts (solid line) and the 
simulated bursts with different detection thresholds for a z peo j. = 2 Rowan-Robinson star forming rate. In particular, the dotted and dashed 
lines are simulated distributions with fiuence threshold of 5xl0~ 8 and lxl0~ 5 erg cm -2 , respectively, (b) Same with (a), but for z pea k = 1. 
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Fig. 2. — (a) Redshift vs. flucnce plot for the simulated bursts (dots) and the observed bursts (crosses) from (Bloom et al. 2003) with a 
z peak = 2 Rowan-Robinson star forming rate, (b) Same with (a), but for z pea k = 1. 
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Fig. 3. — Normalized distributions of the observed jet angles (solid line) obtained from Bloom ct al. (2003) and the simulated jet angles 
with detection threshold of 5xl0~ 8 (dotted line), 5xl0 -7 (dashed line), and lxlO -5 erg cm -2 (dash-dotted line), respectively. The reason 
to choose a high fluence limit, such as lxl0~ 5 erg cm~ 2 , is that only a small fraction of total bursts have their jet break angles measured. 
The median fluence of the 16 bursts with jet break angle measurements is 2.3xl0 -5 erg cm -2 (Bloom et al. 2003), even higher than the 
highest fluence threshold we adopt. 
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Fig. 4. — Fluence vs. jet break angle plot for the simulated bursts (dots) and the observed bursts (crosses) from (Bloom et al. 2003) with 
a z peak = 2 Rowan-Robinson star forming rate. 
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Fig. 5. — Redshift vs. jet angle plot for the observed burst sample (crosses) obtained from Bloom et al. (2003) and the simulated bursts 
with detection threshold of 5xl0 _s (dots), 5xl0 -7 (diamonds), and lxlO -5 erg cm -2 (squares), respectively. 
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Fig. 6. — The luminosity function of GRBs with a quasi-universal Gaussian-like jet structure. The simulated luminosity function can be 
characterized by a broken power law, with the power law indices of ~ —2 in the high luminosity part (L > 10 52 erg s — 1 ) and <~ —1 in the 
low luminosity part (L < 10 52 erg s _1 ). This is consistent with that obtained from Schmidt (2001), except for the low luminosity end below 
10 50 erg s _1 . 
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Fig. 7. — log-/V(> P) — logP plots for the BATSE bursts (solid histogram, from Kommers et al. 2000) and the simulated bursts from 
the quasi-universal Gaussian like model (dotted line). The univeral jet fails to reproduce the logA r (> P) — logP distribution (Guetta et al. 
2004a) as it predicts almost a power law log 7V(> P) — log P distribution that unavoidably over-predicts bursts at the low flux end. The 
simulation results from the quasi-univeral Gaussian jet fits well with the BATSE log N(> P) — logP distribution. 



